Soft dipolar spin ice physics and the ordered phase of the frustrated Tb 2 Sn 2 07 pyrochlore magnet 
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From a microscopic model for the pyrochlore antiferromagnet Tb2Sri2 07, including the crystal field Hamil- 
tonian and interactions between the angular momenta, we compute an effective pseudospin-1/2 Hamiltonian 
H e f! that incorporates perturbatively in the effective interactions the effect of excited crystal field levels. We 
obtain the semiclassical ground states of H c fs and find a region of parameter space with a two-in/two-out spin 
ice configuration on each tetrahedron with ordering wavevector q = and with spins canted away from the 
local Ising axes as found in Tb2Sn207. This ground state can also be obtained from a dipolar spin ice model 
in which the Ising constraint is softened. Monte Carlo simulations on the latter model reveal a region of the 
phase diagram with spin ice-like freezing and another with a transition into Tb2Sn207-type long range order. 
We comment on the differences between Tb2Sn2C>7 and the perplexing spin liquid Tb2Ti207. 



PACS numbers: 75.10.Dg 75.10.Jm 75.40.Cx 75.40.Gb 

In recent years, much effort has been devoted to the ex- 
ploration of geometrically frustrated magnetic systems 1]}|4]. 
Among the many such materials, the spin ices H02T12O7 and 
Dy2Ti207 [3] are remarkable for at least two reasons. Firstly, 
they exhibit a low temperature strongly correlated state with 
residual magnetic entropy and no long range order. Secondly, 
the mechanism leading to effectively geometrically frustrated 
interactions in spin ices and the equilibrium low-temperature 
properties that follow from them have been definitively un- 
derstood OUT). Progress on spin ices has been made quanti- 
tatively through the dipolar spin ice model (DSIM) J3j[5j an d, 
conceptually, via an understanding of the self-screening of the 
dipolar interactions 17| and a description of the low -energy 
excitations as deconfined magnetic charges HO. 

That progress on spin ices has been made so steadily can be 
traced partly to the fact that their magnetic moments can be 
treated as classical and Ising-like. This is due to the fortuitous 
smallness of the inter-ion interactions compared to the single 
ion crystal field anisotropy gap. But this is also good fortune 
for those interested in the problem of quantum fluctuations 
in strongly correlated geometrically frustrated systems for, in 
spin ices, one can explore quantum effects by perturbing away 
from the Ising limit by reducing the anisotropy gap. More- 
over, this is not merely an academic problem for it directly 
pertains to the Tb 2 Sn 2 7 (TSO) and Tb 2 Ti 2 7 (TTO) IfTOl 
compounds which both exhibit rich, complex and poorly un- 
derstood behaviors. 

TSO has a 0.87 K transition to a long range ordered state [9 
1 1 1 where the magnetic moments obey the same (two-in/two- 
out) "ice rules" that hold in spin ice materials |3|. However, 
the ordered moments in TSO are slightly canted away from 
the Ising "in/out" directions while displaying fluctuations that 
persist down to the lowest temperature considered. TTO, in 
contrast, despite a Curie- Weiss temperature #cw ^ — 14 K set 
by the magnetic interactions, fails to develop long range order 
down to at least 50 mK 1 10 1, making it a rare example of a spin 
liquid in three dimensions [2]. Despite numerous experiments 



aimed at exposing the essential physics at play in TSO lf9HTTT- 
14 1 and TTO IfTOl fTsHTTj . no microcopic theory has yet been 
able to explain the behavior of these two materials. 

Just as perturbative calculations that treat the de Broglie 
wavelength as a small parameter allow one to describe quan- 
tum corrections to the properties of simple liquids (e.g. ar- 
gon), in the same spirit, we examine in this paper the effects 
of weak quantum mechanical corrections to the DSIM — a 
classical spin liquid of sorts [2j. As a key milestone, we 
consider TSO. By displaying long range order, TSO is par- 
ticularly amenable to conventional experimental probes of its 
physics. This allows us to benchmark our calculations against 
known experimental results. 

Starting from a microscopic model for TSO, we derive a 
low energy effective Hamiltonian, H c g, which incorporates 
virtual crystal field excitations (VCFEs) about the DSIM ifTBl . 
We compute the semiclassical ground states of H c s and find 
three phases in the vicinity of the Ising limit — one of which 
has the same magnetic structure as found in TSO 0. We 
therefore provide a realistic microscopic explanation for the 
magnetic order of TSO. To explore the extent to which TSO- 
like order is generic among pyrochlore oxides with competing 
exchange, dipoles and single ion anisotropy, we consider a toy 
model, H m , with explicit tunable single ion anisotropy. While 
we find that H m admits a TSO-like phase over a wide range 
of parameters, we observe that a spin ice-like state with dy- 
namically inhibited LRO persists close to the Ising limit. By 
exploring the phase diagram over the entire parameter space, 
we identify a fanning out of several different phases upon tun- 
ing away from the classical Heisenberg antiferromagnet spin 
liquid [ 19 1. Finally, on the basis of this work, we offer a novel 
perspective on how TTO relates to TSO. 

Effective Hamiltonian — The simplest microscopic Hamil- 
tonian consistent with the properties of Tb 3+ ions in TSO 
and with an antiferromagnetic 8qw is H = H c { + V where 
H c [ is the crystal field Hamiltonian with parameters taken 
from Ref. [12]. The interaction term, V — H cx + H^d, is 
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the sum of the nearest neighbor isotropic exchange H cx = 
JemYluj) ' between angular momenta J; (J = 6) on 
site i and dipole-dipole interactions Hdd = ^ r nn^2i>j[^i ' 
Jj — 3(Ji ■ Rij)(Jj • Rij)]|Rij|~ 3 . The dipole coupling 
V = (i (g[i B ) 2 /4Trrl n = 0.029 K where the Lande factor 
of Tb 3+ is g = 3/2 and the exchange has been estimated to 
be J ex ~ 0.08 K fUl. R (J = R; - R, where R t is the 
position of the magnetic ion on site i. The nearest neighbor 
distance r nn = 3.69A = a\/2/A [9| and a is the edge length 
of the cubic unit cell. Throughout this paper, the long range 
dipolar interactions are handled via the Ewald method |20|. 

Diagonalizing H c f gives a spectrum of 2J + 1 = 13 states 
with an Ising-like ground state doublet separated from the first 
excited doublet by a gap A ~ 13.8 K fl2l . Because the in- 
teractions, V, are much smaller than A, we treat V as a per- 
turbation. Using degenerate perturbation theory, we derive an 
effective Hamiltonian, H e g, acting in the space 971 spanned by 
the ground doublet on each site and which includes, to leading 
order in V 2 / A, the effect of admixing of the excited crystal 
field levels into the low energy space (i.e. VCFEs) |fl~8 1. Be- 
cause the low energy space is two-dimensional on each site, 
H e ff can be rendered in the form of a pseudospin-1/2 (S e g) 
Hamiltonian d. In the limit V/A -> 0, H cS ~ O(F) is the 
projection of V onto the ground doublet which gives the dipo- 
lar spin ice Ising model (DSIM) 151 [1811211 . Terms to order 
V 2 / A deform H c g away from the Ising limit. 

H e g ground states — To examine the effects of excited 
crystal field levels on the spin correlations, we determine the 
semiclassical ground states of H c g on a cubic unit cell with 
periodic boundary conditions by replacing the pseudospins, 
S c s (Rj), by classical spins of fixed length 1/2. This is akin to 
finding a semiclassical Neel ordered phase as the leading de- 
scription of the broken symmetry phase in a quantum Heisen- 
berg antiferromagnet. The ground state phase diagram as the 
anisotropy gap A and T>/ J cx are varied is shown in panel (a) 
of Fig.[T] For A = 13.8 K corresponding to TSO, the ground 
state for weakly antiferromagnetic exchange J cx < 0.050 K 
(D/Jex > 0.58) has a q = 001 ordering wavevector with 
zero bulk magnetization and the two-in/two-out spin ice rule 
satisfied by the local [111] (Ising) components of S e s on 
each tetrahedron. Three body interactions in H g cause the 
S' c ff(Rj) to cant away from the local Ising directions. The 
"uncanted" variant of this q = 001 long range ordered spin 
ice (LRSIooi) state is one of the two ground states of the 
DSIM for 1/A = [5] (see Fig. [TJ. For J cx > 0.050 K 
(P/Jbx. < 0.58) and A = 13.8 K, instead of the all-in/all-out 
state observed in the DSIM Q (which appears here only for 
1/A < 0.04 K _1 ), we find that the ground state is an ordered 
ice state with ordering wavevector q = and with the ice 
rule satisfied on each tetrahedron (LRSIooo)- Quantum fluc- 
tuations away from these classical ground states are currently 
being investigated and will be reported elsewhere. The all- 
in/all-out to LRSIooo phase boundary, computed for fixed T>, 
has a maximum for small T>/ J ex which indicates that terms 
to order V 2 / A in perturbation theory are not adequate to ac- 



count for the phase boundary for too small A. 

The ground states spin configurations determined from H c g 
are the semiclassical expectation values of the pseudospins, 
S ff. However, the physical observables are, rather, the ex- 
pectation values of the angular momenta (J") for a = x, y, z. 
We therefore compute the (J Q ) from the S" s expectation val- 
ues. We do this perturbatively in the interactions V finding 

(J a ) = (VTV) + (WQ3°"P) + (VrQVV) + ... (1) 

for the expectation values of (J™) within the low energy sub- 
space where Q = J2\i>)ewi 1^) {^\/{ E g ~ E \i>)) and E g is ^ 
ground state energy of the unperturbed ground crystal field 
doublet. The tilde over the components of (J Q ) indicates that 
the components are taken in the local x — y — z frame on each 
sublattice l20l . The expansion in operators on the right hand 
side of Eq. ([T]i can be rendered in the form of pseudospin-1/2 
operators with coefficients that are computed numerically. 

We are mainly interested in the q = ordered ice state 
(LRSIooo)- We compute H c g on a cubic unit cell and consider 
a state vector, within the pseudospin 971 space, corresponding 
to the LRSIooo with canted S s g (R4) that minimizes the classi- 
cal energy. The state vector corresponding to such semiclassi- 
cal ground state is a direct product of canted pseudospins- 1/2 
states on each site of the cubic cell. The expectation value of 
J z , (J z ), is computed from Eq. ([!} on each of the 16 sites as a 
function of J ex , As one would expect from Eq. ([T} the magni- 
tude of (J 2 ) changes as J cx varies. However, the relative signs 
of y on each of the sites are preserved, so the identification 
of the q = ordered spin ice state is unaltered by the operator 
correction of Eq. ([l|. Next, we find the (3 X ) and expecta- 
tion values on each lattice site for which the only contributions 
come from the terms linear in V in Eq. qTj). The result of this 
calculation is that the (J Q ) ground state configuration is the 
same magnetic structure of canted magnetic moments as ob- 
served in TSO O — a structure that we henceforth refer to as 
LRSItso Note that it is crucial that the ordered moments be 
described via the observables (J Q ) in Eq. (fil, as sole consider- 
ation of the expectation values of the S c s(Ri) does not even 
give qualitatively the correct canting direction and disagrees 
with the canting experimentally observed in TSO |9|. 

Soft dipolar spin ice model ~ The essential physics of 
.ffeff is that VCFEs lead to (i) a softening of the otherwise 
Ising spins and (ii) a new (LRSIooo) ordered spin ice state not 
present in the DSIM [5 1. To explore the effects of the softened 
Ising constraint and the extent to which the appearance of a 
LRSIooo state is generic in dipolar pyrochlore systems with 
softened [111] Ising anisotropy, we investigate a toy model, 
H m , with an explicitly anisotropic g-tensor. 

We take H m = H cx + H^d as used in the microscopic 
model H discussed above where the J; are taken to be Heisen- 
berg spins. For this model, we denote the exchange and 
dipole couplings by J and D respectively. We incorporate the 
anisotropy by introducing spins-1/2 S with components de- 
fined via (J s , 3$, J 2 ~) {g^S^g^y^wS 2 ) with 5y > g x 
to achieve an Ising-like anisotropy. Whereas, in H c g, the 
spin softness is generated through VCFEs, in H m we include 
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FIG. 1: (color online), (a) Semiclassical ground state of H e a on 
a single cubic unit cell as a function of 1/A and J^/Ti. The star 
indicates the estimated microscopic parameters of TSO. (b) Ground 
states of the classical spin model H m as a function of couplings D/J 
and the anisotropy X. Ground states are shown for antiferromag- 
netic J and an Ising anisotropy. Note the similar topology of the 
phase diagrams in (a) and (b) close to D/J, T> / ' ~ 0.2 and small 
anisotropy (small 1/A,X). The inset shows the ground states for 
ferromagnetic J for < X < 1. The vertical dashed line indicates 
the line along which Monte Carlo simulations have been carried out. 
The lower panels show the LRSItso ( c ) and LRSIooi (d) spin con- 
figurations on a single cubic unit cell. The vertical (red) arrow in (d) 
shows the [0, 0, 1] direction of the bulk net moment. 



the effects of this physics at the outset. We parameterize the 
anisotropy g = (g±, g±, g\\) by the parameter X such that 
g = i](X, X, 1 — X) with i] chosen so that g is a unit vector. 
When X = 0, (i.e. g± = 0), we obtain a model that cou- 
ples only the Ising components of the S vectors. This model 
for J — was studied in Ref. [22] and a related model with 
D = 0, X = 1/2 and ferromagnetic J, but with a single ion 
anisotropy — T(Sf ) 2 , has been investigated |23|. 

Panel (b) of Figure [TJ displays the semiclassical ground 
states of H m as X and D/J are varied. The ground states 
are computed on a single cubic unit cell using zero tempera- 
ture Monte Carlo. The main figure shows the ground states for 
< X < 0.5 (Ising anisotropy) and antiferromagnetic J > 0. 
The point (D/J = 0,X = 0.5) corresponds to the classical 
Heisenberg antiferromagnet spin liquid on the pyrochlore lat- 
tice fl9l . from which four long range ordered phases fan out. 
For small D/J, the ground state is an all-in/all-out (AIAO) 
state with ordering wavevector q = with no canting of the 
spins away from the Ising directions. For large D/J, one finds 
LRSIooi states with ordering wavevector q = 001 with the 
ice rule satisfied on each tetrahedron. The AIAO state and 
the LRSIooi state are both ground states of the DSIM and the 
boundary between them for X = is consistent with earlier 
findings 0. We find that the canting angle away from the 



Ising directions varies smoothly away from (for X = 0) 
in the LRSIooi state as the parameters vary. In the isotropic 
(X = 1/2, Heisenberg) limit, we recover the Palmer-Chalker 
ground state ll24l . We find that this state remains the ground 
state even away from the isotropic limit with spins lying in 
the local xy planes. At the center of the figure is a dome in 
which the ground state is a LRSItso state with D/J and X 
dependent canting angles. The inset shows the ground states 
for ferromagnetic exchange (J < 0). The LRSIooi ground 
state (for X = 0) extends to finite X (dashed region), where- 
upon it gives way to a q = two-in/two-out ordered state 
that we label LRSIfm- The distinction between LRSItso and 
LRSIfm is that, in the former (latter), the net moment on each 
tetrahedron is smaller (larger) than the moment 4/i/v3 of a 
q = two-in/two-out state with moments /.i constrained to lie 
along the (111) Ising directions. The labelling of the regions 
in the figure is independent of the canting angle away from the 
(111) directions. The canting angle (not indicated) increases 
as X increases with discontinuities at the phase boundaries. 
Within the LRSItso phase, the general trend in the canting 
angle is for it to decrease as D/J increases. 

We have explored the finite temperature phase diagram 
along the line D/J = 0.3 using Monte Carlo (MC) simu- 
lations with parallel tempering. We chose this line because 
it cuts through three sets of ground states as X is varied (see 
panel (b) of Fig.[TJ). For Ising interactions, X = 0, upon low- 
ering the temperature, the conventional paramagnet freezes 
into a spin ice state characterized by a two-in/two-out con- 
straint on each tetrahedron but without conventional LRO. 
While one would expect some residual entropy, in this case, it 
does not coincide with the Pauling entropy of spin ice because 
the spins are classical 3-component spins with a magnetic spe- 
cific heat per spin of ks at T = (see lower right panel of Fig. 
2). Of particular interest, this freezing is observed at finite 
X up to the boundary between LRSIooi and LRSItso- The 
freezing temperatures are indicated in the top panel of Fig. [2] 
for X < 0.07. For X > 0.07 and D/J — 0.3, the simulations 
equilibrate easily with strong evidence of phase transitions 
into long-range ordered phases corresponding to the ground 
states of Fig. [TJ The transition temperatures from the heat ca- 
pacity peak are plotted in the top panel of Fig. [2] The lower 
panels of Fig. [2] shows for X = 0.25, within the LRSItso 
dome, the order parameter F q= o (an order parameter for the 
q = ordered ice [25 1) and the specific heat for different 
system sizes measured in L — the length of the cubic simu- 
lation cell in units of the cubic unit cell edge. Since the sim- 
ulations reveal a freezing transition for X < 0.07 and LRO 
at finite temperature in the (reentrant) LRSIooi region around 
X = 0.4 (see top panel of Fig.|2|l, there exists a crossover be- 
tween a freezing transition (X < X*) and a phase transition 
to LRO (X > X*) within the LRSI 00 i region of Fig.[TJ). We 
have run a series of simulations for L = 2, 3 and 4 finding 
that this boundary lies near X* ss 0.35, largely independently 
of D/ J as far as we can tell. We always observe the onset of 
finite temperature LRO within the LRSItso dome in Fig. [TJd. 

Materials context — Given the considerable microscopic 
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FIG. 2: (color online). Selected Monte Carlo simulation results on 
H-ny. Simulations have been carried out along the line D/J — 0.3. 
The upper panel shows the freezing temperatures (for X < 0.07) 
and transition temperatures {X > 0.07) determined from the heat 
capacity peak for system size L = 3. The ground states as a function 
of X are indicated. The lower left hand panel shows, for X — 0.25, 
the increase of the q = ordered ice order parameter y q =o at the 
temperature is lowered for L — 2, 3, 4. The lower right hand panel 
shows the heat capacity for X = 0.25 and L — 2, 3, 4. 



states) [9|. We used this result to motivate the study of a 
classical spin model H m referred to as a soft dipolar spin ice 
model (SDSIM). This model also exhibits a finite temperature 
transition into LRSItso states over a large portion of its phase 
diagram indicating that it captures the essential physics ob- 
tained from the microscopic model. Therefore, we are confi- 
dent that we have correctly identified the predominant physics 
leading to the q = spin ice ordered phase of Tb2Sn20y and 
the concurrent spin canting away from the Ising directions. 
A corollary of our results is that spin-lattice couplings are 
not necessary to account for the long range ordered phase of 
Tb2Sn2C>7. More generally, the SDSIM reveals the compet- 
ing and important effects of single ion anisotropy, exchange 
and dipolar interactions in rare earth pyrochlore oxides. We 
hope that our toy model H m will serve as a bridge between 
these materials and ultimately help to unravel the fascinating 
phenomena at play in the TL^Ti^Oy spin liquid material. 
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Canada and the Canada Research Chair program (M. G., Tier 
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similarities between TSO and TTO, one wonders whether H m 
might help shed some light on the low temperature collective 
paramagnetism of TTO. Since the anisotropy gap in TTO is 
larger than in TSO, one might consider TTO to have the larger 
Ising anisotropy (smaller X) and to lie in the frozen canted 
spin ice regime of H m . Interestingly, TTO does exhibit hys- 
teresis and a frequency dependent a.c. susceptibility lfl6l[T7l 
below about 350 mK which is consistent with a slowing down 
of the dynamics. However, this glassiness has been suggested 
to come from only a fraction of the moments 1 16 1, the rest ex- 
hibiting fast spin dynamics down to the lowest observed tem- 
peratures consistent with the lack of LRO iflOl . On the other 
hand, low temperature (T < 500 mK) spin dynamics does 
not distinguish the two materials for they are also observed 
in TSO in spite of the LRO [13]. More critically, we note 
that the mechanism underlying the persistent spin dynamics in 
these materials as in other pyrochlores with LRO and gapped 
collective excitations (e.g. Gd2Sn20y [26]) remains an im- 
portant open problem ||4][26). Since H m captures the LRO in 
TSO which remains dynamical, we cannot rule out the possi- 
bility that dynamics may also occur in TTO at T < 350 mK 
alongside spin ice-like freezing as we find in H m . Further 
investigations will be necessary to settle this matter. 

Conclusion — We have shown that a low energy effective 
Hamiltonian H c s for TL^S^Oz derived from a microscopic 
model including crystal field, antiferromagnetic exchange and 
dipolar interactions exhibits semiclassical ground states coin- 
ciding with the magnetic structure of TT^Sr^Oy (LRSItso 
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